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Abstract: In this article a new numerical scheme for dynamical models with applications to chaotic 
systems was developed, we present the analysis of errors in general. The novel numerical method was 
used to address both linear and nonlinear fractional dynamical systems. The technique was used two 
(linear and nonlinear) systems to be solved of the fractional ordinary differential equations as well as 
fours nonlinear chaotic models. 
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1. Introduction 

Differential equations play both active roles in the description of natural and complicated events in 
fluid mechanics, biology, and applied physics, fractional equations it is a powerful tool for describing 
living phenomena in physics and engineering [1-11]. Therefore, numerical scheme has been that has 
gained acceptance as a useful system. Therefore, numerical scheme has been recognized as an effective 
tool for solving systems. We dedicate this article the creation of a novel numerical approach that 
incorporates the fundamental theorem of fractions and the two-stage Lagrangian polynomial [12]. The 
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numerical approach will be used to resolve both linear and nonlinear systems of ordinary fractional 
equations. 

Several methodologies, including physics and engineering, have been applied to problems in 
management, economics and biology, among others [13-18], recent attention has been drawn to the 
fractional SIR Model due to the spread of diseases such as Covid-19. Furthermore, we will apply this 
technique to several chaotic model [19-21], The numerical simulation results shown in both situations 
comprise and three-dimensional phase pictures with varying parameters. The novelty of this study is that 
it provides a numerical solution for fractional derivative order in linear and nonlinear fractional dynamical 
systems. 

The goal of this research was to lay the groundwork for using fractional systems in other fields of 
science and technology. This method's value resides in the fact that it can be applied to a wide variety of 
models from disease models and chaos models to additional models in pathology and dynamical models 
in order to locate a numerical solution. 

2. Preliminaries 

Definition 1 Riemann-Liouville fractional integral operator of order @ > 0 for a function y(z) is 
given by [22]: 


1 t eee oe 
D* yt) = a [ @ —a* * 1y"@dz = I* *y"@,t > 0 (1) 
Definition 2 The fractional integral of order @ of a function f is defined as [23] 
I; y(t) eat y(t) + = J, y@dt.t=>00<a<1 (2) 


Definition 3 The Mittag-Leffler function is a generalization of the exponential function. This 
function can be expressed as follows: 


E,(t) = Uke=0 FaKs4) (3) 


Definition 4 For yeH1(0,T), T > 0,@e(0,1] Then the AB fractional operator [24] y(t) in the 


Riemann-—Liouville is given by 
ABN ee — @_+z¢ 
Dry® = — J, y@Ea ‘= (t—1) ) dr. O<a<1 (4) 
In this expression B(@) satisfies the condition B(0) = BQ) = 1. 
Definition 5 For y € H! (0, T), T > OThen the AB fractional operator [24] y(t) in the Caputo sense is 


given by 


AB _ Bla) ptay() a 

Diy) = —— J, Fa (Gas (t- D*) dt.0<a<1 (5) 
In this expression B(@) satisfies the condition B(0) = B(1) = 1. 
3. Applications: 


3.1 The new numerical scheme 
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© F(t,.x(t,)) 


eo a— 
Xne1 = Xp + ABC(a) 


= h* f(t Xx) ee | | . 
aad (“Ee T(a+ 2) ((@—k+ 1"(n—k+2+a)—(n—b)"(n—k + 2+ 2a) 


_ hf (ty4.X4-1) 
ne + my 


BS ey (t,)) 


((n— k+ 1)" —(n—k)*(n—k4+14+ «)) 


Vane Yyt oe ) 


a h°f(t..y,) > | P | 
aad. <? T(a + 2) ((n—k+1)*(n—k+ 2+ a@)—(n—k)*(n—k+ 24 2a)) 


WF (ta.%_4) 


T(a +2) (@—k+ 1" —(@—k*(n-k +14 ) 


®) F(t, ,2(t,)) 


Zn41 = Ygt eae ) 


a Rh “f(te zy) - - 
aa °F ea T(@+2) ((n4+1—k)*(n—k+2+ a) —(n—k)*(n—k +2 + 2a)) 
h* f (ty_1,Z4_1) +1 
— —k+1)*"—(m—kK*(m—k4+14+ 
r+ 2) ((n—k+1) (n—k)*(n—k+1+2)) 
The use of the novel numerical approach [16] for solving fractional differential equations is covered 
in the next section. 


3.2. linear systems of the fractional ordinary differential equations: 

DAS = aS, — bSz (6) 

Das = cS, — dS, 

We apple AB operators in system (6-7) we get: 

ABDAS = aS, — bSz 

ABDAS = cS, — dS; (7) 

Where a=2, b=1, c=4, d=3. 

We simulate the numerical solution using the specified numerical approach., the estimation of error 


is given 
Table 1: Estimation of Error 
siiedei < Rall = lanes ~ Vexac(t) ! 
a,= 0.94 10 0.4875 
a = 0.97 100 0.033 
a;=0.99 1000 0.0057918680 
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Fig.2 Comparison of numerical solution and exact solution for h = 0.01. 
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Fig.3 Comparison of numerical solution and exact solution for h = 0.001. 


3.3. Nonlinear systems of the fractional ordinary differential equations: 
S; = O(S2 + Sy) 

S2 = —Sz — 5153+ BS, 

S3 = —S3+ SS; (10) 

We apple AB operators in system (6-7) we get: 

APD?S, = a(S, 5 i S,) 

AB D?S2 = —Sz — S1S3+ BS; (11) 

ABD}S3 = —S3 + S452 


Where 0 <A=1. 
Equilibria: The three equilibrium points obtained — by 


Ey — (Js - LJP = 1,8 — 1), E, = (0,0,0), E = (-yeF-L—~jp-if 1) 


Jacobean matrix of the system (8) is given 


—6 7] 0 
J(E,) = G =a) 2 -s,| (12) 
Ss S,; -1 


At the parameter values @ = 14,8 = 110,the equilibrium points are 
E, = (10.45, 10.45, 109),(13) FE; = (10.45,10.45 ,109), E, = (0,0,0) 
The Jacobean matrix at (10.45 , 10.45 , 109) 


solving (8) are 


—15 15 0 
I(E,) = 1 —1 -1045 (14) 
10.45 10.45 =i 
The characteristic equation obtained from (6) at equilibrium E; ( S; ,55 ,S3 ) is as in the following 
equation. 
w* + 17? + 125 w+(3270=0 (14) 


The solution of the characteristic equation gives the eigenvalues. similarly, the eigenvalues 


corresponding to E, and E3 are calculated. 
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Table 2: Equilibria of (11) and their eigenvalues 


E; Equilibria Eigenvalues Stability 
E, (10.45 , 10.45 , 109) (1.15074+12.96i, - 19.31) — Unstable 
E, (0,0,0) (33.2189, -1, -49.22) Unstable 
E; (—10.45 ,—10.45 , 109) (1.1507+412.96i, -19.31) Unstable 


Fig.4. Distribution equilibrium points of the fractional-order system. 
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Fig.5. Stability of the fractional-order system. 
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Fig.6. Chaotic attractor for a = 0.98, h=0.land t=2000 in a y-x plane. 


Fig.7. Chaotic attractor for a = 0.98, h=0.land t=1500 in a y-x plane. 


Fig.8. Chaotic attractor for o = 0.98, h=0.land t=500 in a y-x plane. 
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Fig.9. Chaotic attractor for o = 0.98, h=0.land t=500 in a y-x-z space. 


4. Conclusion: 

We proposed a new numerical approach to solve fractional differential equations (Linear and 
nonlinear) made up from this kind of derivative. The novel scheme technique combines polynomials and 
Lagrange theorems with fractional calculus' basic theorem. When compared to precise answers, the 
suggested technique can be demonstrated to be extremely accurate. The method is simple, effective and 
can be widely applied as it shows chaos, the stability analysis of this method has discussed with error 
estimation. Effective and adaptable to many fractional systems, this technique is widely used. 
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